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The present paper concerns the derivation of phase-integral quantization conditions for the two- 
centre Coulomb problem under the assumption that the two Coulomb centres are fixed. With this 
restriction we treat the general two-centre Coulomb problem according to the phase-integral method, 
in which one uses an a priori unspecified base function. We consider base functions containing three 
unspecified parameters C, C and A. When the absolute value of the magnetic quantum number m 
is not too small, it is most appropriate to choose A = \m\ =^ 0. When, on the other hand, \m\ is 
sufficiently small, it is most appropriate to choose A = 0. Arbitrary-order phase-integral quantiza- 
tion conditions are obtained for these choices of A. The parameters C and C are determined from 
the requirement that the results of the first and the third order of the phase-integral approximation 
^ \ coincide, which makes the first-order approximation as good as possible. 

In order to make the paper to some extent self-contained, a short review of the phase-integral 
method is given in the Appendix. 

. PACS numbers: PACS numbers: 03.65.Sq, 31.15.-p, 31.15.Gy 

in 

pH . 1. INTRODUCTION 

The two-centre Coulomb problem, that is, the problem of solving the Schrbdinger equation for the motion of an 
electron with the charge — e (e > 0) in the field of two fixed Coulomb centres with charges Z\e and Z^e at the 
distance r\i from each other, plays an important role in several fields of theoretical physics, for example in the theory 
^ , of diatomic molecules, in scattering theory, and in the three-body problem. The two-centre Coulomb problem has 
therefore been the subject of extensive studies both by numerical and by approximate analytical methods, and hence 
the literature on this problem is very comprehensive. In spite of this fact it is still of interest to continue the treatment 
of this problem for arbitrary values of Z\,Zi and r\i- One reason for this is that the numerically exact solution of 
the problem meets with difficulties when \Z\ — increases. There appear also numerical difficulties for large values 
of r\2- For a general review of the problem we refer to Eyring et al. Herzberg M, Slater || and Rosen H. 

Ignoring the finiteness of the mass of the protons, Bates et al. [pP obtained important numerical results for the 
""^5 I hydrogen molecule ion. Corresponding numerical results were obtained by Wallis and Hulburt || for the homonuclear 
. — ■ one-electron two-centre problem, by Bates and Carson H for the ion HeH 2+ [Z\ = I, Zi = 2), by Wind Q for 
the ground state of the hydrogen molecule ion, by Peek || for the states lsa g and 2pa u of the hydrogen molecule 
ion, and by Ponomarev and Puzynina jl~C, O for several states of the Coulomb two-centre system with Z± = 1 
and Z2 — 2,3, ...,8. Hunter and Pritchard|f 2j used a numerical procedure to compute nonadiabatic energies for the 



> 



first few rotation- vibration levels of S s i?2 , YiHD + and HgD^ ■ For the hydrogen molecule ion Rosenthal and E. 
Bright Wilson, Jr., pj| calculated an accurate value of the ground state energy which is in agreement with the values 
obtained by Wind Q and Peek||. For different internuclear distances Bates and Reid Murai |T^] and Murai and 
Takatsu |ll| calculated electronic energies for various states of the hydrogen molecule ion under the assumption of 
infinite proton mass. For a large range of internuclear distances Winter et a/.p"7|] made very accurate calculations for 
the lowest 20 states of the molecule ion HeH 2+ . Klaus [p"8| studied the electronic energy of the ground state of the 
hydrogen molecule ion for small internuclear separation r\i and confirmed the remarkable discovery by Byers Brown 
and Steiner Jl9| that the electronic energy cannot be expanded in powers of r\2 alone, but that powers of lnri2 must 
also be included. Klaus jl8| also obtained for the hydrogen molecule ion further terms in the series given by Byers 
Brown and Steiner [ ft9| for the general two-centre Coulomb problem. 

There exist also various approximate analytical methods for solving the quantal two-centre Coulomb problem. In 
one of them one uses the quasiclassical approximation, that is, the first order of the phase-integral approximation. In 
using this method, one encountered in the early papers difficulties associated with the divergence of the phase-integral 
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due to an inappropriate choice of the phase- integrand (also called the quasimomentum) , which is not determined 
quite uniquely. Different authors have mastered these difficulties in different ways, but no single unifying method has 
until now been proposed for the two-centre Coulomb problem. An essential feature in our method is the use of the 
phase-integral approximation generated from an unspecified base function; see the Appendix. 

The semiclassical quantization of the low-lying electronic states of the hydrogen molecule ion was treated by Strand 
and Reinhardt (20). Pajunen J21J calculated the energy levels of the hydrogen molecule ion (under the assumption 
of infinitely heavy nuclei) in the first and the third order of the phase-integral approximation. Although he refers to 
one of the papers (his reference 10) in which the phase-integral approximation generated from an unspecified base 
function was introduced, he does not use the freedom to choose optimal expressions for the functions that he denotes 
by Qmod(n) and Q m od( v ), an d that in the present paper correspond to the more general base functions Q(£) and Q(jf), 
respectively. In the present paper we shall make full use of the possibility to choose Q(^) and Q(rj) most conveniently. 

The phase-integral method, in which one uses the phase-integral approximation generated from an unspecified base 
function [ p2[ , offers a method for mastering the previously mentioned difficulties in a unified way for an arbitrary 
order of the phase-integral approximation. In the present paper we shall apply this method to the quantal two-centre 
Coulomb problem with fixed Coulomb centres. 

For the convenience of the reader a short review of the phase-integral method and formulas to be used are given in 
the Appendix. 

The present paper will be the basis for further work, in which convenient transformations to complete elliptic 

and Lakshmanan et aZ.[^4|, p5[ . 
Thereby the contour integrals, occurring in the quantization conditions, will be expressed in terms of complete elliptic 
integrals. 



integrals will be used, analogously as in papers by Lakshmanan and Kaliappan|23 



2. SEPARATION OF THE SCHRODINGER EQUATION IN ELLIPTIC COORDINATES 

We start by quoting some well-known results. The time-independent Schrodinger equation for the motion of an 
electron of mass [i and charge — e (e > 0) in the field of two fixed Coulomb centres with charges Z\e and Z%e takes 
the following form 

^A,-^-^W)=W), (2.1) 
2^ ri r 2 J 

where r\ and r 2 are the distances of the electron from the two centres, r is the position vector of the electron, and E is 
the electronic energy. To obtain the total energy one must add the potential energy of the two fixed charges, getting 

Etotal = — — — — h E, (2.2) 
r\2 



where ri2 is the distance between the two centres. The differential equation (2.1) is separable in elliptic coordinates. 
If one introduces the variables 

n + r-} 

£=— -, 1<£<+oo, (2.3a) 

r\2 

V= H^I1 7 -l<rj<+l, (2.3b) 

and puts 

*(r) = X(e)r(ry) e ™^, (2.4) 

where m is the magnetic quantum number (positive or negative integer or zero), and <f> is the corresponding angle, 
the separation yields, in atomic units (h = e = fj, = 1), the two differential equations 

^ ' 2\ ^ \ i I 1 2 



*, v (1 - ) + r ? + bv - A - w ' Y = 0> (2 - 5b) 
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where A is the separation constant and 



Putting 



P 2 = -\ r l2E, (2.6a) 

b' = r 12 (Z 2 + Z x \ (2.6b) 
b = r 12 (Z 2 -Z 1 ). (2.6c) 



*(0 = 7t/%t, (2.7a) 



Y (v) = (2.7b) 



we can transform the differential equations (|2.5|) into 



^+R(o)f(0 = o, (2- 



flfa) fffa) = 0, (2.8b) 



where 



br] — A' m 2 — 1 
1 -t? 2 ~ (1 - rff 



, in - 4' m 2 - 1 , 
R(V) = -P 2 + -i T - T, ™ > ( 2 - 9b ) 



with 

A' = A-p 2 . (2.10) 



The differential equations (2.S) are of the Schrodinger type. By solving them simultaneously under the boundary 
conditions that /(+1) = /(+00) = and g(— 1) = g(+l) = 0, one can obtain the energy and the separation constant 
as functions of the distance r% 2 and the quantum numbers. 

When Z\ = Z 2 every bound-state wave function Y(rj) is either an even or an odd function of 77, and when r is 
reflected at the centre of symmetry for the two-centre Coulomb proble m, th e wave function ( |2.4[ ) remains unchanged 
when Y(rf) and m are both even or both odd, while the wave function (2.4) changes sign when one of Y(r)) and m is 
even and the other is odd. 



3. APPLICATION OF THE PHASE-INTEGRAL METHOD 



The essential features of the phase-integral method are briefly sketched in the Appendix. The phase-integral 
solutions of the differential equations (2.8a) and (2.8b), respectively, are linear combinations of the phase-integral 
functions 

r*(0«p|±</%(0^| (3-la) 

and 

g _ *(?7)exp j±« J q{rj)dr\, (3.1b) 

respectively, where q(£) and q(rf), respectively, are generated according to (A5a,b), (A6a,b,c), (A3) and (A2) in 
the Appendix, with R(z),Q(z) replaced by the appropriate functions i?(£),Q(£) and R{f]), Q(r]), respectively, the 
quantities pertaining to the ^-equation being characterized by a tilde. 
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3.1. Base functions generating phase-integral solutions 

As is seen from (2.9a,b), the functions i?(£) and R{rj) have poles at £ = ±1 and rj = ±1, respectively; these poles 
are of the second order if to ^ ±1, but of the first order if m = ±1. Furthermore, we note that when m ^ the 
coefficients of the second-order poles differ from j, while for to = they are equal to \. 

There are two main alternatives, discussed in the Appendix for the case of the radial Schrodinger equation, for 
choosing the base functions generating the phase-integral functions (3.1a) and (3.1b). Unifying these two alternatives, 
we choose the squares of the base functions to be 

*a/« 2 A' — C + b'£ A 2 . . 

Q 2 (0 = -P 2 + g2 _ 1 - ^3^, (3.2a) 



\ 2 | -A' + C + b v A 

W = -P + : 2 n 2W ( 3 - 2b ) 

1 — r/ z (1 — r) ) 



where C, C and A are parameters, the choice of which we shall discuss below. The introduction of these parameters 
increases essentially the flexibility of the phase-integral formulas to be obtained. 

By choosing C = C = 1/4 one obtains in the limit when m — > the energy E and the reduced separation constant 
A' correctly from the first-order phase-integral quantization conditions (to be derived in Section 3.2). For arbitrary 
values of r\2 it is most appropriate to determine C and C as functions of ria such that one obtains the same value 
of p 2 (i.e., of the energy) and also of A' in the first and in the third order of the phase-integral approximation. One 
thereby achieves the optimal accuracy obtainable from the first-order quantization conditions. We emphasize that 
this can be achieved by the use of the phase-integral approximation generated from an unspecified base function 
(described in the Appendix), but cannot be achieved by means of the JWKB approximation. The decisive properties 
of the phase-integral approximation in question versus the JWK B ap proximation have been explained in some detail 
by Dammert and P.O. Froman |^6f and by Froman and Froman p2f . 

When \m\ is not too small, we choose A = \m\ ^ 0, but when \m\ is sufficiently small we choose A = 0. For to = 
one should always choose A = 0. The two alternatives A = \m\ ^ and A = yield solutions with different properties. 
When A = \m\ ^ the phase integral solutions of the differential equations (2.8a) and (2.8b) remain valid as £ — > ±1 
and i] — > ±1, respectively. When A = the phase- integral solutions of the differential equations in question break 
down as £ — > ±1 and rj — > ±1, respectively, but the regular solutions of the £- and ry-equations can be obtained at 
some distances from those points by the use of the connection formula described in subsection A.2.b of the Appendix, 
when one there replaces I by (|to| — l)/2; cf. (2.9a,b). The wave functions obtained in that way are expected to be 
the more accurate the stronger the Coulomb singularities of R(£) and R(r/) at £ = ±1 and rj = ±1 are. However, even 
if the Coulomb singularities are strong, these wave functions are not expected to be good if \m\ is too large, in which 
case one should use A = \m\ ^ 0, as mentioned in the beginning of this paragraph. 

Using the terminology classically allowed region and classically forbidden region in a generalized sense, viz. to 
characterize regions where Q 2 (£) or Q 2 (ri) is larger than zero and less than zero, respectively, we shall now discuss 
the wave functions pertaining to the ^-equation and the ^-equation. 



3.1.1. Wave functions pertaining to the ^-equation 

According to (3.2a) the function — Q 2 (£), where 1 < £ < +oo, always corresponds to a single-well potential. 

For A = |m| ^ there are four zeros of Q 2 (£), which we denote by £1,^2, £3 and £4; see Fig. 1(a). The zeros £1 
and £2 may be real and both less than 1, or they may be complex conjugate. The zeros £}(> 1) and £3) are 

real. There is thus a classically allowed region between £3 and £4, but classically forbidden regions for 1 < £ < £3 and 
£ > £4. The phase-integral wave functions generated from Q(C) are good at £ = 1, and we can use the arbitrary-order 
connection formula, given by (A13) and (A14) in the Appendix, for tracing the physically acceptable wave function 
from the classically forbidden region between 1 and £3 to the classically allowed region, that is, the region between £3 
and £ 4 (> £3) where Q 2 {£) > 0; see Fig. 1(a). 

For A = there are two zeros of Q 2 (£,), which we denote by £s(< £4) and 1). The classically allowed region 
lies between £3 and £4 when £3 > 1, but between 1 and £4 when £3 < 1; see Fig. 1. The phase-integral wave function 
is not good at £ = 1, but for the physically acceptable wave function one can obtain a phase-integral expression in 
the interior of the classically allowed region, when £3 > 1 and the classically forbidden region between 1 and £3 is 
sufficiently large, by means of the arbitrary-order connection formula given by (A13) and (A14) in the Appendix, and 
when £3 < 1 by means of the connection formula presented in subsection A.2.b of the Appendix. 
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Both when A = \m\ ^ and when A = the wave function, obtained in the classically allowed region to the left of 
£4 as described above, can be joined to the physically acceptable wave function traced from the classically forbidden 
region to the right of £4 into the classically allowed region to the left of £4 with the aid of the arbitrary-order connection 
formula given by (A13) and (A14) in the Appendix. In this way alternative quantization conditions, corresponding to 
A = \m\ ^ and A = 0, can be obtained. They can be combined into one quantization condition. Fig. 1 illustrates 
the two possible situations that the classically allowed region is delimited either by £3 and £4 [Fig. 1(a)] or by the 
pole at £ — +1 and the turning point at £ = £4 [Fig. 1(b)]. 



3.1.2. Wave functions pertaining to the r)-equation 

The function —Q 2 (rj) may correspond either to a single- well potential or to a double- well potential. 

For A = I to I 7^ the phase- integral solution is valid at the poles r\ = ±1 (which delimit classically forbidden regions) 
and can be traced into the classically allowed region closest to the pole in question with the aid of the connection 
formula given by (A13) and (A14) in the Appendix. When there is only one classically allowed region, one obtains 
the quantization condition by identifying the two expressions for the wave function in that region. This case applies 
when —Q 2 (j]) is a single- well potential (Fig. 2), or when the energy lies so far above the top of an underdense barrier 
[Fig. 4(b)] that it is appropriate to disregard the presence of the complex conjugate zeros 7/2 and % of Q 2 (rj). When 
—Q 2 {v) corresponds to a double- well potential, the wave function can be traced from the region on one side of the 
barrier to the region on the other side with the aid of the arbitrary-order connection formula for a barrier described in 
Section A. 3 of the Appendix; see Fig. 3 and Fig. 4. Joining the two expressions for the wave function thus obtained 
to each other, one obtains the quantization condition. 

When A = [Fig. 2 or Fig. 4(b)] the phase-integral wave function is not good at 77 = ±1, but at some distance 
from these points physically acceptable solutions can be obtained in the same classically allowed region by the use of 
the connection formulas presented in Sections A. 2 and A. 3 of the Appendix. Thus one obtains two expressions for 
the wave function, and by identifying these expressions one obtains a quantization condition. 



3.2. Quantization conditions 

3.2.1. Quantization conditions pertaining to the ^-equation 

For the differential equation (2.8a) the physically relevant interval is 1 < £ < 00. The phase-integral quantization 
condition for the situation in Fig. 1(a) involves a contour integral in the complex £-plane encircling £ 3 and £ 4 , while 
the phase- integral quantization condition for the situation in Fig. 1(b) involves a contour integral encircling the 
simple pole at £ = 1 and the generalized classical turning point £4. The quantization condition (A31) in the Appendix 
applies to the first situation, and the quantization condition (A32) in the Appendix applies to the second situation. 
Introducing the notations 

N 



L = Y J L (2n+1 \ (3.3a) 

n=0 

£(2„+l) 1 f -(2n+l) ( £ )c (33b) 



2 



N 



L' = J2L {2n+1 \ (3.4a) 



n=0 



L^ n +^ = - [ <7 (2n+1) (£R, (3.4b) 

where q(£) is obtained according to (A5a,b), (A6a,b,c), (A3) and (A2) in the Appendix, and A^ and A^, are the 
appropriate contours of integration pertaining to £3 > 1 and £3 < 1, respectively, and shown in Fig. 1(a) and Fig. 
1(b), respectively, we can write the two quantization conditions corresponding to £3 > 1 and £3 < 1 as follows 

.5+-Jtt, &>1, (3.5a) 
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where s is an integer. If, when A = \m\ ^ and hence £3 > 1, we enlarge the contour of integration in Fig. 1(a), 
so that the new contour A^, encloses the turning points £3 and £4 as well as the pole at £ = 1, and if we compensate 
in (3.5a) along with (3.3a,b) for this change by taking the residue of the integrand at £ = 1 into account, we obtain a 
general quantization condition, valid for both cases £3 > 1 and £3 < 1, i.e., for both situations depicted in Fig. 1, viz. 

L'= (jT+ s +l)*> £ 3 >lor6<l. (3.6) 

Besides condensing the two alternative quantization conditions (3.5a) and (3.5b) nicely into one formula, the quanti- 
zation condition (3.6) has the further merit that, if the integration along the contour is made numerically (in cases 
where expressions in terms of complete elliptic integrals are not available), it may be advantageous to use the contour 
A 1, instead of A^ when £3 > 1. 

The quantization condition (3.6) yields the value of the reduced separation constant A' as a function of p 2 and C; 
see (3.2a). 



3.2.2. Quantization conditions pertaining to the r\-equation 

In the physically relevant interval — 1 < 77 < 1 the function —Q 2 (r\) may correspond to a single- well potential (Fig. 
2) or to a double-well potential with a superdense (Fig. 3) or underdense (Fig. 4) barrier. When, in the case of an 
underdense barrier, the energy lies sufficiently far above the top of the barrier, it may be preferable to disregard the 
barrier and to treat the double-well potential problem as a single-well potential problem. 

When A = \m\ 7^ or A = 0, and ~Q 2 {f]) is or can be considered as a single- well potential, and the classically 
allowed region is delimited by two simple zeros of Q 2 (rj), as shown in Fig. 2 and in Fig. 4(a), we obtain from (A31) 
the single-well quantization condition 



L = 



■S+-I7T, 



s = non-negative integer, 



(3.7) 



where by definition 



2 



q(r])dr), 



(3.8) 



At, being a closed contour encircling the generalized classical turning points. Note that in the derivation of (3.7) we 
have considered the classically forbidd en r egions to be thic k also when A = (Fig. 2). When A = |m| ^ we can 
with the aid of residue calculus write (3.7) along with (3.8) as 



L'= (| m |+ s +~)7r, (3.9) 

where L' is defined by 

L ' = \l l(v)dv, (3.10) 

J A L , 

A li being a closed contour encircling —1 and +1; see Fig. 2 and Fig. 4(a). 

When A = 0, and the classically allowed region is delimited by two first-order poles of Q 2 {rf), as shown in Fig. 
4(b), and the energy lies far above the top of the barrier, one can consider —Q 2 {r)) as a single-well potential. From 
(A33) one then obtains the quantization condition (3.9) with L' defined by (3.10), where Aj,> is now the contour in 
Fig. 4(b). 

We disregard the possibility that A = and the residues of Q 2 {rf) at 77 = —1 and 77 = +1 have different signs, since 
this case has so far not appeared in the applications. 
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When —Q 2 (rj) corresponds to a double-well potential (Fig. 3 and Fig. 4), which is usually the case, the quantization 
condition (A39) gives 

cos(a + + 4> 2a) = ^-l) (3 . n) 
[1 + exp(— 27rK)\ 2 



where 

n. = 

2 



a = £ when A = \m\ j= 0, (3.12a) 



o=(|m| + l)^ when A = 0, (3.12b) 



N 

a = Y / a (2n+1) , (3.13a) 

n=0 

a (2n+l) = Re l f g (2n+l)^ d ^ ) (3 13b ) 

2 J A„ 



AT 



/3 = ^/3 (2n+1) , (3.14a) 

n=0 

/3 (2„+l) = _ Rc l j q (2n+l)^ dT] ^ (3Mb) 
2 As 



N 



K = Y,K 2n , (3.15a) 



n=0 



K2n = ^-[ q {2n+1 Hv)dv- (3.15b) 
For the super-barrier case (Fig. 4) we can instead of (3.15b) use the alternative formula 

K 2n = -2Im±-J^ q^ n+1 \ri)dr, 

= _2Im— / q^+^^drt, (3.156') 

which is useful in connection with the transformation to complete elliptic integrals. With the use of (3.13a,b), (3.14a,b) 
and (3. 15a, b') each one of the quantities a, (3 and K is then obtained as the real or imaginary part of an integral 
over the contour A Q or kp. The contours of integration A Q , A^ and Kk for sub-barrier and super-barrier energies are 
shown in Fig. 3 and Fig. 4. The analytic expressions for the quantity <j> are given in terms of K and K 2n by (A28) 
and (A29a,b,c) in the Appendix. The quantity (f> is of decisive importance for energies in the neighbourhood of the 
top of the barrier. The quantities a and /3 are positive. The quantity K is positive when 7i 2 and 773 are real, it is 
equal to zero when 772 and 773 coincide, and it is negative when 772 and 773 are complex conjugate. The quantity K 2 
may be positive or negative irrespectively of whether 772 and 773 are real or complex conjugate. 

Analogously as we changed the original contour A^ into A^,, when dealing with the quantization condition (3.5a), 
we can, when A = |m| 7^ 0, change the contours A Q and A^ depicted in Fig. 3(a) and Fig. 4(a), so that, instead 
of letting them enclose only 771 , 772 and 773, 774, respectively, we make each one of them enclose also a pole, 77 = —1 or 
fj = +1, respectively Calling the new contours k a > and Ap> [see Fig. 3(a) and Fig. 4(a)] and defining 

N 

a' = J2^' {2n+1 \ (3.16a) 

n=0 
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a '(2n+l) = Re l f g (2n+l)^ dj7j (3 lgb) 

2 "'A„, 



AT 



/3' = ^/3' (2 " +1) , (3.17a) 



n=0 



/3 '(2„+l) = _ Rc i /" g (2n+l)^ d)7) (317b) 

2 

and recalling that the functions Y 2n are regular analytic at r\ = ±1 when A = |m| 7^ 0, we find with the use of residue 
calculus that for A = |m| 7^ 

a' - a = - /3 

= Att/2, (3.18a) 

a' -13' = a - P 

With the aid of (3.18a,b) we obtain from (3.11) along with (3.12a) for A = \m\ 7^ the following quantization condition 

cos[«' + (3' + 4>-{\m\ + 1).] = - ^^L , (3.19) 

[1 + exp(— 27ta)J 2 

which has the same merits, relative to the original form of the quantization condition, i.e., (3.11) with (3.12a), as were 
mentioned in connection with the quantization condition (3.6), pertaining to the ^-equation. When A = we define 

N 

a' = Y,a {2n+1 \ (3.20a) 

n=0 

a '(2n+l) =a (2n+l) j (3 20b) 



13' = Y J P (2n+1 \ (3.21a) 



n=0 



/3 '(2n+l) =/3 (2n+l) ) (321b) 

and note that (3.18a) is obviously valid also for A = 0. With the use of the theory of complex integration one finds 
that also (3.18b) is valid for A = 0. Then it follows from (3.11), (3.12b) and (3.18a,b) that the quantization condition 
(3.19) is valid also when A = 0. The quantization condition (3.19) thus covers in a unified and convenient form both 
cases A = \m\ 7^ and A = 0. 

The quantization condition (3.19) can be rewritten as 

a' + 0' + cj>- (|m| + 1)tt 

cos[67r/ (2p)l „ , 

= iarccos l - — 1 v y 'l , + 2s'tt, 3.22 

[l + exp(-2 7 r J ft')]2 



where s' is an integer. As already mentioned, (3.18b) is valid for A = |m| 7^ as well as for A = 0. When the plus sign 
in (3.22) applies, we use (3.18b) to express (3' in terms of a', and when the minus sign in (3.22) applies, we use (3.18b) 
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to express a' in terms of (3'. Replacing s' by s a or sp, we thus obtain from (3.22) and (3.18b) the two quantization 
conditions 

ml 1\ 6 bn 

1 cos (%) 

+ -arccos — — — r , (3.23a) 

2 [l + cxp(-27rX)]5' 

1 cos (%) 

— arccos v ' _ . , (3.23b) 

2 [l + cxp(-2^)]5' 

where we choose the branch of arccos such that the last two terms on the right-hand side of (3.23a,b) cancel in the 
limit K — ► +oo, that is, when the barrier becomes infinitely thick. In this limit, the formulas (3.23a,b) simplify to 

a '=(jy + s « + I) ^ ( 3 - 24a ) 
P=(jr + 'f> + \)«- ( 3 - 24b ) 

For the particular case that we have a symmetric two-centre Coulomb problem, i.e., that Z\ — Z 2l as is the case for 
the ion , the double- well potential pertaining to the ^-equation becomes symmetric (6 = 0), and the quantization 
conditions (3.23a,b) can be simplified: 

H \\ 4> 

-T + S « + 2) 7T -2 
+ iarctan exp(— irK), (3.25a) 

P = {— + S " + 2) 7r -2 

— iarctan exp(— nK). (3.25b) 

The reduced separation constant A' , obtained from (3.23a,b) in the general case and from (3.25a,b) in the symmetric 
case, is a function of p 2 and C. 



3.2.3. Comments on the quantization conditions 

In the existing semi-classical treatments, the quantization conditions derived on assumptions valid for m ^ are 
in general extrapolated to m = (corresponding to a particularization of our case A = 0) without any motivation. 
It is, however, not allowed to obtain a quantization condition corresponding tom = from a quantization condition 
corresponding to m ^ by letting m take continuous values and tend to zero by a limiting procedure. That the 
formulas obtained by such an extrapolation are valid is an a posteriori conclusion. The correct justification of the 
quantization conditions for A = rests on the use of the connection formula in subsection A.2.b of the Appendix. 

For given values of ri2, m and E one can, as already mentioned, obtain the possible values of the reduced separation 
constant A' in the differential equation (2.8a) with (2.9a) by applying phase-integral quantization conditions for a 
single-well potential, while to obtain the possible values of the reduced separation constant A' in the differential 
equation (2.8b) with (2.9b) one has to use quantization conditions either for a single-well potential or for a double- 
well potential. The appropriate quantization condition for the ^-equation determines A' as a function of p 2 and C. 
The appropriate quantization condition for the ^-equation determines A' as a function of p 2 and C. The eigenvalues 
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of p 2 , and hence the energy eigenvalues E, are obtained from the requirement that these two expressions for A' must 
be equal to each other. One then obtains A' from the quantization condition for the ^-equation or the ry-equation. 
The value thus obtained for A 1 depends obviously on the choice of C and C. One should choose these parameters 
as funcitons of m in such a way that very accurate values of p 2 and A' are obtained already in the first order of 
the phase-integral approximation. A practically useful criterion for this can be formulated as follows. For every 
value of 7*12 one determines C and C such that the first-order approximation gives the same value as the third-order 
approximation for p 2 (and hence for the energy) as well as for A' . In this way one can make all calculations within 
the framework of the phase-integral method. 
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APPENDIX A: PHASE-INTEGRAL METHOD 

Since the present paper is based on phase-integral formulas that are scattered in different publications, we collect 
in this Appendix the background material that is necessary for reading the paper. 
The phase-integral method for solving differential equations of the type 

+ iZ(z)V- = O (Al) 

involves the following items: 

1. Arbitrary-order phase-integral approximation generated from an unspecified base function Q{z), as described in 
Chapter 1 of p2] ]; see also Dammert and P.O. Froman pc| . 

2. The method for solving connection problems developed by Froman and Froman |2jj |, generalized to apply to the 
phase-integral approximation referred to in the above item. 

3. Supplementary quantities, expressed analytically in terms of phase-integrals. An example is the quantity </>, 
which is a new notation for the quantity — 2a in Froman et a/. |28f| , and which is of decisive importance, when 
two transition zeros lie close to each other. 

We shall first briefly describe the phase- integral approximation referred to in item 1. Then we collect connection 
formulas pertaining to a single transition point [first-order zero or first-order pole of Q 2 (z)] and to a real potential 
barrier, which can be derived by means of the method mentioned in item 2 combined with comparision equation 
technique for obtaining the supplementary quantity <fi mentioned in item 3 and appearing in the connection formula 
for a real barrier. Finally we present quantization conditions for single- well and double- well potentials, which can be 
derived by means of the connection formulas just mentioned. These quantization conditions are used in our treatment 
of the two-centre Coulomb problem. 

1. Phase-integral approximation generated from an unspecified base function 

For a detailed description of the phase-integral approximation generated from an unspecified base function we refer 
to chapter 1 in [ p2| . A brief description is given below. 

In the arbitrary-order phase-integral approximation in question there appears an unspecified function Q(z) called 
the base function. This function is often chosen to be equal to Ri(z), but in many physical problems it is important 
to use the possibility of choosing Q(z) differently when there exist certain exceptional points, for example the origin 
in connection with the radial Schrodinger equation, and, correspondingly, the poles of Q 2 (£) and Q 2 {r]) at £ = 1 and 
77 = ±1 in the two-centre Coulomb problem. In the present paper we introduce in the base functions a parameter A, 
chosen such that either A = \m\ 7^ or A = 0, and two parameters C and C to be determined such that the first- and 
third-order results coincide, in order that the first-order approximation be as good as possible. 
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To be able to write the phase-integral approximation generated from an unspecified base function in condensed 
form one introduces the new independent variable 



C= / Q(z)dz (A2) 

and the function 



_ i , , R(z) — Q 2 (z) 
e = Q-Hz)^Q-H Z ) + LL_|li . (A3 ) 



It can be shown that in a local region of the complex z-plane, where the absolute value of £o is small, the differential 
equation (Al) has the approximate solutions 

ip = q~^(z)exp[±iw(z)], (A4a) 



w(z) = / q(z)dz, (A4b) 

where the lower limit of integration zq is an unspecified constant, and the function q(z), pertaining to the phase-integral 
approximation of the order 2N + 1, is given by 

N 



q{z) = Y,1 [2n+1) ( Z )' ( A5a ) 
n=0 

q^+^iz) = Q(z)Y 2n , (A5b) 

with the first few functions Y 2n given by 

Y = l, (A6a) 

Y 2 = ~ £o , (A6b) 

*~H-±£- (A6c) 

The choice of the function Q(z) does not affect the expressions for Y<z n in terms of £q and £; only the expressions for 
£o and C as functions of z are affected. 

It is an essential advantage of the phase-integral approximation described above versus the Carlini|2^] (JWKB) 
approximation in higher order that the former approximation contains the unspecified base function Q(z), which one 
can take advantage of in several ways. A criterion for the determination of the base function is that the function 
So be in some sense small in the region of the complex z-plane that is relevant for the problem under consideration. 
However, this criterion does not determine the base function Q{z) uniquely; it turns out that, within certain limits, 
the results are not very sensitive to the choice of Q{z) 1 when the approximation is used in higher orders. With a 
convenient choice of Q(z) already the first-order approximation can be very good. On the other hand an inconvenient, 
but possible, choice of Q(z) introduces in the first-order approximation an unnecessarily large error that, however, in 
general becomes corrected already in the third-order approximation. 

The freedom that one has in the choice of the base function Q{z) will be illuminated in a concrete way below. For 
a radial Schrodinger equation the usual choice of Q 2 (z) is 

Q 2 (z)=R{z)-^. (A7a) 



However, the replacement of (A7a) by 



2 / „ I const 



Q\z) = R(z) - - —, (A7b) 
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where the coefficient of \ should be comparatively small, does not destroy the great accuracy of the results usually 
obtained with the phase- integral approximation in higher orders. There is thus a whole set of base functions that 
may be used, and there are various ways in which one can take advantage of this nonuniqueness to make the choice 
of the base function well adapted to the particular problem under consideration. For instance, by adapting the choice 
of Q 2 (z) to the analytical form of R(z) one can sometimes achieve the result that the integrals occurring in the 
phase-integral approximation can be evaluated analytically. To give an example we assume that R(z) contains only 
exp(z) but not z itself. In this case it is convenient to replace the choice (A7b) by the choice 

^o/ s . . 1 const 

Q 2 W = *M ~ i(-F3Tp - — V ( A7c ) 

By a convenient choice, for instance of the unspecified coefficient in (A7b) or (A7c), one can sometimes attain the result 
that, for example, eigenvalues or phase-shifts are obtained exactly for some particular parameter value in every order 
of the phase-integral approximation. By making this exactness fulfilled in the limit of a parameter value, for which 
the phase-integral result without this adaptation would not be good, one can actually extend the region of validity of 
the phase-integral treatment; see p. 12 in |^0| . When the differential equation contains one or more parameters, the 
accurate calculation of the wave function may require different choices of the base function Q{z) for different ranges 
of the parameter values. To illustrate this fact we consider the radial Schrodinger equation. For sufficiently large 
values of the angular momentum quantum number I we obtain an accurate phase-integral approximation (valid also 
close to z = 0) if we choose Q 2 (z) according to (A7a) or (A7b). If the value of / is too small, this phase-integral 
approximation is not good. It can be considerably improved (except close to z = 0), when the absolute value of the 
coefficient of \ in R(z) is sufficiently large, if one chooses instead 

Q 2 ( z ) = R( z )+ l l±±l. (A7d) 

The corresponding phase-integral approximation is not valid close to z — 0, but the wave function that is regular and 
tends to z l+1 as z — > can be obtained sufficiently far away from z = by means of the connection formula that will 
be presented in subsection A.2.b of this Appendix. 

The appearance of the unspecified base function Q(z) in the phase-integral approximation is thus very important 
from several points of view. In our treatment of the two-centre Coulomb problem we use two essentially different 
kinds of base function [corresponding to A = \m\ ^ and A = in (3.2a,b)], which yield approximate solutions with 
different regions of validity. 

When the first-order approximation is used, it is often convenient to choose the constant lower limit of integration 
z in (A4b) to be a zero or a first-order pole of Q 2 (z). This is, however, in general not possible when a higher-order 
approximation is used, since the integral in (A4b) would then in general be divergent. If z Q is an odd-order zero or 
an odd-order pole of Q 2 {z), it is therefore convenient to replace the definition (A4b) of w{z) by the definition 



w(z) = - / q{z)dz, (A8) 

where T za (z) is a path of integration that starts at the point corresponding to z on a Riemann sheet adjacent to 
the complex z-plane under consideration, encircles z a in the positive or in the negative sense, and ends at z. It is 
immaterial for the value of the integral in (A8) if the path of integration encircles zq in the positive or in the negative 
sense, but the endpoint must be z. For the first-order approximation the definitions (A4b) and (A8) are identical. 
It is useful to introduce a short-hand notation for the integral in the right-hand member of (A8) by the definition 

/ q(z)dz = \ f q(z)dz. (A9) 
J(z ) l Jr H) (z) 

For the first order of the phase-integral approximation one can replace (zq) by zq in the left-hand member of (A9) 
and thus get an ordinary integral from zq to z instead of half of the integral along the contour T Zo (z). In analogy to 
(A9) one defines a short-hand notation for an integral in which the upper limit of integration is an odd-order zero or 
an odd-order pole of Q 2 (z). When one has two transition points of that kind as limits of integration, the definition of 
the short-hand notation with both limits within parentheses implies that the integral is equal to half of the integral 
along a closed loop enclosing both transition points. The simplified notation in the left-hand member of (A9) for the 
integral in the right-hand member of (A9) was introduced by Froman et al. pl| , pp 160-161. It makes it possible 
to use, for an arbitrary order of the phase-integral approximation, a similar simple notation and almost the same 
simple language (although in a generalized sense) as for the first order of the phase-integral approximation. One thus 
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achieves a great formal and practical simplification in the treatment of concrete problems, when an arbitrary order 
of the phase-integral approximation is used. 

We remark that the notations used above differ from the notations in the original papers published up to the 
beginning of the eighties in the respect that Q 2 {z) and Qmod( z ) m th° se papers correspond in the present paper to 
R(z) and Q 2 (z), respectively. 



2. Connection formulas associated with a single transition point 



a. Connection formulas pertaining to a first-order, real transition zero 



Before the phase-integral approximation generated from an unspecified base function had been introduced, N. 
Froman ]3^j derived arbitrary-order connection formulas associated with a turning point for the particular phase- 
integral approximation of arbitrary order corresponding to Q 2 (z) — R(z). After the phase-integral approximation 
generated from an unspecified base function had been introduced, it turned out that these connection formulas remain 
valid also when Q 2 (z) ^ R(z). Below we shall present the general connection formulas. 

The functions R{z) and Q 2 (z) are assumed to be real on the real z-axis (the a;-axis). On this axis there is a 
generalized classical turning point t, i.e., a simple zero of Q 2 (z), and there is, in a generalized sense, a classically 
allowed region on one side of t, i.e., a region where Q 2 {x) > 0, and a classically forbidden region on the other side of 
t, i.e., a region where Q 2 (x) < 0. Defining 

w{x) = [ q(z)dz, (A10) 

•A*) 

we can write the connection formula for tracing a phase-integral solution of (Al) from the classically allowed to the 
classically forbidden region as 



.4 



q~i(x) expji \w(x) \ + ^ } 

q~*{x) expj-i \w{x)\ + ^ j — 
— >(A + B) q~^(x) exp[\w(x)\] , 



-B 



B 



where A and B are constants, which are arbitrary except for the requirement that 
zero. As a consequence of (All) we have the connection formula 



cos 



sin 5 



\w(x)\ 
q~?(x) exp[|w(a;)| 



(All) 

must not be too close to 



(A12) 



where S is a real phase constant that must not be too close to zero. The connection formula for tracing a phase-integral 
solution of (Al) from the classically forbidden to the classically allowed region is 



(x) exp[—\w(x)\]+C q 2 (x) exp[|w;(a;)| 



q *(x) 



cos 



|u>(x) 



and it is valid provided that the condition 



Cexp [|w(a:)|] < exp [— |i/>(x)| 



(A13) 



(A14) 



is fulfilled. For a numerical studyof the accuracy and the properties of the connection formula (A13) with C = we 



refer to N. Froman and Mrazek|33 . We emphasize the one-directional character of the connection formulas (All), 
(A12) and (A13), which means that the tracing of a solution must always be made in the direction of the arrow. This 
pro perty of the connection formulas has been thoroughly investigated and even illustrated numerically by N. Froman 
|]34j for the first order of the Carlini HjJ (JWKB) approximation. The whole discussion in [H applies in principle 
to the connection formulas for the higher orders of the phase-integral approximation as well. The above connection 
formulas for the phase-integral approximation of any order may in many cases be used for obtaining very accurate 
solutions of physical problems, when the classical turning points are well separated, and when there are no other 
transition points near the real axis in the region of the complex z-plane of interest. Within their range of applicability 
the connection formulas are very useful because of their simplicity and the great ease with which they can be used. 
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b. Connection formula pertaining to a first-order transition pole 



Now we assume that in a certain region of the complex z-plane around a first-order transition pole t, i.e., a first-order 
pole of Q 2 (z), we have 

R(z) = — 7^— — rl H h a regular function of z, (A15) 

[z — ty z — t 



D 



t 



a regular function of z, 



(A16) 



where 21 + 1 is a non-negative integer. We assume that the absolute values of B and B are sufficiently large, while 
the absolute value of B — B and the absolute value of the difference between the two regular functions in (A15) and 
(A16) are sufficiently small. There is one particular curve on which w(z), defined as 



w(z) 



(0 



q{z)dz, 



(A17) 



is real. For the first order of the phase- integral approximation this is the anti-Stokes line that emerges from t. Therefore 
we use, also when a higher order of the phase-integral approximation is used, the terminology an anti-Stokes line that 
emerges from t in a generalized sense to denote the anti-Stokes line on which w(z), defined in (A17), is real. For the 
first-order approximation (and under certain unnecessarily restrictive assumptions) Froman and Froman p7| obtained 
a phase-integral formula [their eq.(7.28)], valid sufficiently far away from t on the anti-Stokes line that emerges from 
t, for the particular solution ip(z) of (Al) that fulfils the condition 



lim . 

z^t [z - t) 



1+1 



1. 



(A18) 



This formula can be generalized to be valid for an arbitrary order of the phase-integral approximation generated from 
an unspecified base function and can then be formulated as follows. On the lip of the anti-Stokes line emerging from 
t, where arg w(z) is smallest, the solution of (Al) that fulfils the condition (A18) is, sufficiently far away from t, given 
by the phase-integral formula 



irc- 



w{z) 



» (z)cos 



\w(z)\ 



1 + 



(A19) 



where c is the residue of 2 at z — t and is thus determined by the expansion of ip(z) in powers of z — t, and 



the sign of 7TC |^j has to be chosen conveniently. For the special case that I = one finds that c = B. Formula 

(A19) can easily be particularized to the case that R(z) and Q 2 (z) are real on the real z-axis (the x-axis) and t lies 
on that axis. 



3. Connection formula for a real, smooth, single-hump potential barrier 

Our starting point is a paper by Froman and Froman [ B5| . Although it was assumed in the treatment there that 
Q 2 (z) = R{z) with the notations in the present paper, the results obtained are valid also when Q 2 (z) ^ R(z). In the 
present paper it is convenient to introduce partly other notations than in pB| . Thus we now denote by t' and t" the 
two relevant zeros of Q 2 (z), i.e., the two generalized classical turning points in the sub-barrier case (t 1 < t") and the 
two complex conjugate transition zeros in the super-barrier case (Im t' < 0,1m t" > 0). Let x' , called x\ i n [R5[ , be 
a point in the classically allowed region of the real z-axis to the left of the barrier, and let x" , called X2 in [p5|, be a 
point in the classically allowed region of the real z-axis to the right of the barrier. We introduce the notations 

6 = \F 22 \ exp(K), (A20a) 



■d = arg F 22 , 



(A20b) 
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= -2a = 



arg F 12 , 



(A20c) 



where F 12 and F 22 are denned in Q, and K is defined by eq.(12) in |3jJ. In the definitions (A20a,b,c) it is assumed 
that the phase of q^{z) is chosen as shown in Fig.l in j35|. We shall, however, in the following write the formulas in 
such a way that they remain unchanged if one changes the phase of (z) and hence the phase of <p (z); see (A5a,b). 
The quantity K in (A20a) is then given by 



K = 



1 



2' 



q(z)dz 

A 

(t") 

1 I q(z)dz, 
(?) 



(A21) 



where A [not to be confused with the parameter A in the base function Q(r))] is a closed contour of integration 
encircling both t' and t" , but no other transition point, and the integration is performed in the direction that in the 
first- order approximation yields K > for energies below the top of the barrier and K < for energies above the top 
of the barrier. If higher-order approximations are used, the quantity K may become negative also for energies below 
(but not too far from) the top of the barrier; see Table I in |3(|. We have replaced a, defined by eq. (28) in |35|, by 

— I [cf. (A20c)] in order to get better agreement with a notation used by other authors; see for instance Child fj7| . 
Now we define 



B' = A : exp 



(A22a) 



A' = B 1 exp(+ij), 



(A22b) 



where the notations in the right-hand members are those used i n J35[ . Using the short-hand notation defined in (A9), 
we obtain from eqs. (25a) and (9a), with x\ replaced by x' , in |35[ and (A22a,b) in the present paper 









r' 


ijj(x') = A' 


q-i(x') 


exp 


Re / q(z)dz 








J(t>) 









r' 




B' 




exp ^— i 


Re / q(z)dz 


) 








J(t>) 





(A23a) 

and from eqs. (25b) and (23), with x 2 replaced by x" , in |35| and (A20a,b,c) and (A22a,b) in the present paper 









r" 


i){x") = A" 


q-i(x") 


exp 


Re / q(z)dz 








J(t») 



where 









r" 




+B" 


q-i(x") 


exp ^— i 


Re / q(z)dz 


) 








J(t") 





B") ~ ^ \B' ' 



(A23b) 



(A24) 



~ ( 0exp [-z(f + d)\ {9 2 + l)5exp( +l 0; 
l2 + l)5exp(- i 0) 6>exp [+i(f +1?)] 



(A25a) 



detM 



(A25b) 



It is seen from (A23a,b) that the coefficients A' and A" are associated with waves travelling away from the barrier, 
while the coefficients B' and B" are associated with waves incoming towards the barrier. From (A25a,b) we obtain 



= ( ^cxp [-i(f - 0)] (6 2 + l)5exp(+^; 
\(9 2 + l)Texp(-i^) Sexp [+t(f - ■&)] 



(A26) 
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One thus obtains ill -1 from M by replacing d by — d. We emphasize that the above formulas are in principle exact, 
provided that one knows the quantities 9, i9 and </>, which depend slightly on x' and x" . Furthermore, the two transition 
zeros associated with the potential barrier need not lie far away from other possibly existing transition points; see 
p5| . However, when one introduces for 6, $ and 4> the approximate values that will be given below, the barrier is 
assumed to be well separated (in the sense just described) from all transition points that are not associated with the 
barrier. 

When A 1 and B' are given constants, associated with a wave function that is given at the point x', the coefficients 
A" and B" , which are obtained from (A24) along with (A25a), depend slightly on x 1 and x" via the quantities 9, $ 
and </>. One obtains the derivatives of ijj(x') and ip(x") from (A23a,b) by considering A',B',A" and B" formally as 
constants. 

When the transition points that are not associated with the barrier lie sufficiently far away from t' and t" , it follows 
from eq.(43a) in H) and (A20a,b) that 

6t*sexp(K), (A27a) 



i?«0. (A27b) 

The quantity <j> is particularly important when the energy is close to the top of the barrier, but it is important also 
for energies well below the top, if one with the use of higher orders of the phase- integral approximation wants to obtain 
very accurate results. In practice one cannot obtain useful expressions for <f> from the exact formula (A20c). Under 

the assumption that d is not too close to zero at the top of the barrier, Froman et al. ]3§| derived by means of 
comparision equation technique, adapted to yield formulas for supplementary quantities in the phase-integral method, 
an approximate, but very accurate, formula in the (2N + l)th order of the phase-integral approximation [their eqs. 
(5.5.30), (5.5.25a-g), (5.4.23) and (5.4.21)], from which we can obtain the formula 



are 



N 



rgr Q + iKj - K ln\K \ + J2 (2 " +1) , (A28) 

where 



n=0 



s(3) 



(1) = K , (A29a) 

= (A29b) 

24^0 



with 



4>™ = * +-5* *L, (A29c) 

2880^ 2AK% 2K K J 



K 2n = -L / Y 2n Q(z)dz, n = 0,1, 2, ....N, (A30a) 
2tti J a 



N K 

i 7T 



(A30b) 



A being the previously described contour of integration encircling t' and t" but no other transition point, with the 
integration performed in the direction that makes Kq positive when t' and t" are real, i.e., when the barrier is 
superdense, but negative when t' and t" are complex conjugate, i.e., when the barrier is underdense. (Note that we 
perform the integrations in (A21) and (A30a) in opposite directions in order to make these formulas to agree with 
eq. (12) in |i5| and eq. (5.4.21) in |||, respectively.) The result given by (A28), (A29a,b,c) and (A30a,b) can also 
be obtained from p8| , where —2a is the same as our (f>. 

We emphasize again that for the validity of (A28) with the expressions (A29a-c) for 0( 2n+1 ) the essential restriction 



is that 



dp 



must not be too small at the top of the barrier, which means that close to its top the barrier is 



approximately parabolic. However, when the energy is close to the top of the barrier, it is the slight deviation from 
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parabolic shape close to the top of the barrier that determines the values of the quantities -ft^n, n > 0, and one 
needs accurate values of these quantities for obtaining accurate values of cf> in higher orders of the phase-integral 
approximation. 

The barrier connection formula presented in this section is valid uniformly for all energies, below and above the 
top of the barrier. We would also like to emphasize that while the connection formulas pertaining to an isolated 
turning point (N. Froman |32|) are one-directional, the barrier connection formula (A23a,b) with (A24) and (A25a), 
which is valid when the barrier is well isolated and has an approximately parabolic top, is bi-directional. However, 
the neighbourhood of an energy that corresponds to a resonance requires a careful discussion. 



4. Quantization conditions for single-well and double-well potentials 

In this section we shall present quantization conditions for general single-well potentials |5^, [l0|, 41 and double- 
well potentials j2f| |4l], 4^, 43 1, valid for any conveniently chosen order of the phase-integral approximation, in forms 
especially adapted to the two-centre Coulomb problem. In the quantization conditions pertaining to the double-well 
potential there appears the supplementary quantity (j>, which was discussed in Section A. 3, and which is of particular 
importance for energy eigenvalues in the neighbourhood of the top of the barrier; cf. numerical results in ^8|, where 

Q 2 , Q 2 mod and a correspond to R, Q 2 and — J, respectively, in the present paper. Comparision with numerical results 
p8[ shows that all energy eigenvalues, also the low-lying ones and those in the neighbourhood of the top of the 
barrier, are obtained very accurately from the phase-integral quantization conditions when the third- or fifth-order 
approximation is used. 

Since arbitrary-order phase-integral quantization conditions for the single-well and for the double-well potential 
problem have been given in previous work, we restrict ourselves to quoting those results and taking into account the 
fact that we are dealing with the special potentials pertaining to the two-centre Coulomb problem. 



a. Quantization conditions for single-well potentials 

We assume that R(z) and Q 2 (z) are real on the real z-axis (the x-axis) and that there are two transition points t' 
and t" (> t 1 ) on this axis, each one of which may be either a first-order transition zero, i.e., a first-order zero of Q 2 (z), 
or a first-order transition pole, i.e., a first-order pole of Q 2 (z). These transition points are assumed to lie far away 
from all other transition points. On the real axis between t' and t" it is assumed that Q 2 (x) is positive. With the aid 
of the connection formulas in Section A. 2 we can derive the quantization conditions that will be presented below. 

When both t' and t" are first-order transition zeros, we obtain the quantization condition p% M 



{t") 



(*') 



q(z)dz 



8+ \ 



7T, S = 0,l,2. 



(A31) 



When one of the transition points t' and t" is a first-order transition zero, and the other is a first-order transition 
pole in the neighbourhood of which R(z) and Q 2 (z) can be expanded according to (A15) and (A16) with / = g -1 , 
i.e., 1(1 + 1) = m t~~— , we obtain under the assumptions introduced below those expansions the quantization condition 



(*") 



q(z)dz 



(//) 



|m| 
2 



S+ 2 



(A32) 



where m is an integer (positive, negative or zero), and s is also an integer. 

When both transition points t' and t" are first-order transition poles, in the neighbourhood of which R(z) and Q 2 (z) 
can be expanded according to (A15) and (A16) with I = ^ m ^~ 1 , i.e., 1(1 + 1) = m ~ 1 for both transition poles but 
possibly with different coefficients B and B for the two transition poles, we obtain under the assumptions introduced 
below (A15) and (A16) the quantization condition [B0| 



(*") 



(*') 



q(z)dz 



1 

S+ 2 



7T, 



(A33) 



where m is an integer (positive, negative or zero), and s is also an integer. 
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b. Quantization conditions for double-well potentials 

We assume that R(z) and Q 2 (z) are real on the real z-axis (the x-axis) and that there are either two generalized 
classical turning points t 1 and t" (> t') [real, simple zeros of Q 2 {z)] associated with a superdense potential barrier or 
two complex conjugate simple transition zeros t' and t" [simple zeros of Q 2 (z); Im t' < 0, Im t" > 0] associated with 
an underdense potential barrier. The classically allowed region to the left of the barrier is to the left delimited by 
a transition point t_ (on the real axis), and the classically allowed region to the right of the barrier is to the right 
delimited by a transition point t + (on the real axis), where t- and t + are both either generalized classical turning 
points, i.e., simple zeros of Q 2 (z), or first-order transition poles, i.e., first-order poles of Q 2 (z). The points i_ and 
t+, as well as other possibly existing transition points that are not associated with the barrier, are assumed to lie far 
away from t' and t" . When t_ and t+ are simple transition poles we assume that in the neighbourhood of t± we have 

R(z) = — H h a regular function of z, (A34) 



'(*) 



B± 

z — t- 



a regular function of z, 



(A35) 



where m is an integer (positive, negative or zero), B± and B± are sufficiently large to their absolute values, while the 
absolute values of B± — B±, as well as the difference between the regular functions in (A34) and (A35), is not too 
large. Under the above assumptions we obtain with the aid of connection formulas in Section A. 2 and Section A. 3 
the quantization condition 



where 



tan ( 


At') 

Re / q(z)dz 
J{t~) 


* \ ( 
+ — — a 1 tan 1 


At") 
Re / q(z)dz 
J(t+) 






[l + exp(-2K)]i 


- 1 




[l + exp(-2i<r)]3 


+ 1' 




K = 


— / q(z)dz, 
2« J A 





(A36) 



(A37) 



A being a closed contour, enclosing t' and t", along which the integration is performed in the direction that makes the 
first-order contribution to K positive when the barrier is super-dense but negative when the barrier is under-dense, 
and 



■| when t_ and t + are transition zeros, 
(|m| + 1)^ when t_ and t + are transition poles. 



(A38) 



The quantization condition (A36) can be rewritten into the form 



Re 



(*') 



(*-) 



q{z)dz 



Re 



(*") 



q(z)dz 



(t+) 



2a 



Re j£* \ q(z)dz — Re jA* ) q(z)dz 



[l + exp(-2 J ftT)]2 



(A39) 
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FIG. 1: The figure gives schematic pictures of — Q 2 (£) for £ > 1 and of the contours of integration in the complex £-plane. 
The cuts are indicated by wavy lines. On the upper lip of the cut Q(£) is positive. Only those zeros of <3 2 (£) that are relevant 
for the contours of integration are shown. For A = \m\ ^ (a) always applies, and the relation between the integrals associated 
with the contours A^ and A^, is L' = L + the zeros £i and £2 of Q 2 (£), which are not shown in the figure, may be real or 
complex conjugate. For A = there are only two zeros, £3 and £4, of Q 2 (£), and either (a) or (b) may apply. 
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FIG. 2: The figure gives schematic pictures of — Q 2 (ri) for —1 < 77 < 1 and of the contours of integration when — Q 2 {rf) in the 
interval — 1 < r\ < 1 is a single-well potential, which may occur for A = \m\ / as well as for A = 0. The cut is indicated by a 
wavy line, on the upper lip of which Q{rf) is positive. The contour A L i can be used only when A = \m\ 7^ 0, and the relation 
between the integrals associated with the contours Al and A L i is then L' — L + \m\. The quantization conditions, expressed 
in terms of complete elliptic integrals, for the situation in this figure with A = \m\ 7^ are the same as the corresponding 
quantization conditions for the situation in Fig. 4(a). The quantization condition for the case in the present figure with A = 
has so far not appeared in the applications. 
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FIG. 3: The figure gives schematic pictures of —Q 2 (r]) for —1 < 77 < 1 and of the contours of integration in the complex 
»)-plane, when —Q 2 (rj) corresponds to a double-well potential with a super-dense barrier. In (a) A = \m\ 7^ and in (b) A = 0. 
The cuts are indicated by wavy lines, and the parts of the contours of integration that lie on Riemann sheets adjacent to the 
complex jy-plane under consideration are dashed. In the left-hand classically allowed region Q(rj) is positive. Only those zeros of 
Q 2 (rj) that are relevant for the contours of integration are shown. We recall the relations (3.18a) which mean that a' = a + 4p 
and/3'=/3+4p. 





FIG. 4: The figure gives schematic pictures of — Q 2 (»?) for —1 < r\ < 1 and of the contours of integration in the complex 7/-plane, 
when —Q 2 {ri) corresponds to a double-well potential with an under-dense barrier. In (a) A = \m\ ^ and in (b) A = 0. The 
cuts are indicated by wavy lines, and the parts of the contours of integration that lie on Riemann sheets adjacent to the complex 
77-plane under consideration are dashed. In the left-hand classically allowed region Q{rf) is positive. Only those zeros of Q 2 (ri) 
that are relevant for the contours of integration are shown. If the energy lies sufficiently far above the top of the barrier, one 
may treat the double-well potential problem as a single-well potential problem with the classically allowed region between 771 
and 7/4 in (a) and between the poles at 77 = — 1 and 77 = +1 in (b). One introduces then new cuts [between r/i and 7/4 in (a) and 
between —1 and +1 in (b)], on the upper lips of which Q(rj) is positive, and uses the contours Al and A L i . In (a) the relation 
between the corresponding integrals is L' = L + \m\. In both (a) and (b) L is related to a and (3 by the relation L = a + (3. 



